Effective charge and free energy of DNA inside an ion channel 
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Translocation of a single stranded DNA (ssDNA) through an a-hemolysin channel in a lipid 
membrane driven by applied transmembrane voltage V was extensively studied recently. While the 
bare charge of the ssDNA piece inside the channel is approximately 12 (in units of electron charge) 
measurements of different effective charges resulted in values between one and two. We explain these 
challenging observations by a large self-energy of a charge in the narrow water filled gap between 
ssDNA and channel walls, related to large difference between dielectric constants of water and lipid, 
and calculate effective charges of ssDNA. We start from the most fundamental stall charge q s , which 
determines the force F s — q s V/L stalling DNA against the voltage V (L is the length of the channel). 
We show that the stall charge q s is proportional to the ion current blocked by DNA, which is small 
due to the self-energy barrier. Large voltage V reduces the capture barrier which DNA molecule 
should overcome in order to enter the channel by |qfc|V, where q c is the effective capture charge. 
We expressed it through the stall charge q s . We also relate the stall charge q a to two other effective 
charges measured for ssDNA with a hairpin in the back end: the charge q u responsible for reduction 
of the barrier for unzipping of the hairpin and the charge q e responsible for DNA escape in the 
direction of hairpin against the voltage. At small V we explain reduction of the capture barrier with 
the salt concentration. 



I. INTRODUCTION 

A DNA molecule in a water solution carries negative 
charges. With the help of applied voltage, it can translo- 
cate through a wide enough ion channel located in a 
lipid membrane [H, 0, 0, [J 0, 0, 0, 0, Hi or through a 
solid state nanopore in a semiconductor film [llj, ITll ll2j . 
An intensively studied example is the translocation of a 
single stranded DNA (ssDNA) molecule through an a- 
hemolysin (a-HL) channel 0, 0, 0, II 0, 0, 0, 0, H ■ With 
the average internal diameter ~ 1.7 nm the channel is 
wide enough for ssDNA molecules, but is too narrow for 
a double helix. To be specific below we always talk about 
the experimental data for this system. Our theory is also 
applicable to a double helix DNA translocating through 
a narrow nanopore [T3, EH EH , but there is less quanti- 
tative data for this case. 

In order to study the translocation experimentally, the 
electric current through the channel is observed under 
voltage V, applied between two vessels of salty water on 
both sides of the membrane. Due to large conductivity 
of the bulk solution practically all the voltage drops on 
the membrane. When a ssDNA is added to the negative 
voltage side it is dragged into the channel by the voltage 
(Fig. [I}. When the ssDNA molecule is in the channel 
as shown in Fig. [T] the ion current is blocked, and the 
blocked current lb is much smaller than the open pore 
current Iq without DNA in it, 
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Translocation events in a single channel can be studied 
by monitoring the current. It was argued recently that 
the steric mechanism of strong current blockage is ampli- 
fied by the increase electrostatic self-energy of an ion in 
water passage narrowed by DNA Q because DNA and 
lipids have much smaller dielectric constants than water 



and, therefore, the electric field lines of an ion in the 
space between DNA and the lipid are squeezed in the 
channel increasing the ion self-energy This idea was bor- 
rowed from the physics of narrow ion channels without 
DNA 0. 

Besides the blocked current, one can measure the time 
between two successive translocation events, r, or the 
capture rate R c = 1/r of DNA molecules into the chan- 
nel. It is natural to compare the observed value of R c 
with the diffusion limited rate Rd of ssDNA capture. 
This comparison shows that R c -C Rd- For instance 0], 
the typical R c is in the range of 0.01 - 10 s- 1 at applied 
voltage 50 — 200 mV and ssDNA concentration 0.9 fiM, 
while Rd ~ 100s _1 . The ratio R c /Rd may be as small as 
10~ 6 if one extrapolates the experimental data to V — 0. 
So there must be a large barrier ~ \AksT for ssDNA cap- 
ture. We return to the nature of this barrier in the end 
of Introduction, but first we concentrate on challenging 
question of the voltage effect on this barrier. 

The capture rate at zero voltage R c (0) is so small that 
all experiments are actually done with a large applied 
voltage V = 50 - 200 mV . The voltage pulls DNA into 
the channel and reduces the barrier for DNA capture. It 
was found [l|, y, |J] that the capture rate is 



Rc(V) = R c {0)ex P (~q c V/k B T), 



(2) 



where q c — — 1.9e is an effective "capture" charge, and 
e is the proton charge. Apparently for ssDNA in a-HL 
channel \q c \ is much smaller than the absolute value of 
the total DNA charge in the channel —N^e ~ — 12e and 
this is why large voltage V ^> ksT/e is necessary Q in 
order to make the capture rate observable. 

Why is the capture charge of DNA so small and what 
does it depend upon? How is the capture charge related 
to the stall charge defining the stall force F s , which one 
should apply to DNA occupying the whole length of the 
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FIG. 1: The side view of the membrane and the channel with 
captured DNA. All DNA phosphates in the shaded part are 
neutralized by K + ions. The contact layer on each end has 
the length D. The arrow shows the direction of DNA motion. 



channel (Fig. [T]) to stall it against the voltage VI (This, 
for example, can be done with the help of an laser tweez- 
ers [HI, EH). I n other words, F s = —F p , where F p is the 
force, with which the voltage V pulls the stalled DNA. 
We write F s as 



-q s E, 



(3) 



where the electric field E = —V/L and q s is the stall 
effective charge. The charge q s seem to be the simplest 
and the most fundamental effective charge one can intro- 
duce for DNA. Is it different from q c ? If yes, which one is 
larger? How are these two charges related to unzipping 
and escape charges which describe ssDNA with a hairpin 
at the end (see definitions below)? 

Inspired by all these challenges in this paper we use 
for DNA the simplest model of a rigid cylinder charged 
by the point like surface charges (phosphates located on 
the spiralling backbone) and moving coaxially through 
a cylindrical tunnel filled by salty water (Fig. [1]). The 
model of rigid cylinder should be good for a double helix 
DNA in a narrow cylindrical semiconductor pore. For 
ssDNA such model gets some support from the known 
tendency of ssDNA to stuck its bases [l5[ in a bulk so- 
lution. It is natural to expect that this tendency is en- 
hanced inside the channel. One may say that the case of 
ssDNA in cv-HL channel pushes our model too close to 
the molecular limit. Nevertheless, we will show that our 
results for effective charges are in a reasonable agreement 
with experiment. We believe that this happens because 
these results are practically model independent. 

The plan of this paper is as follows. In Sec. [IT] we 
discuss our model and the electrostatics of the channel. 
We concentrate on the role of the enhanced self-energy of 
a salt ion in the narrow water-filled space between DNA 
and internal walls. First, we argue that DNA in the 
channel is almost perfectly neutralized and, second, salt 
cations are bound to DNA charges. Then we introduce 
narrow charged contact layers near the end of DNA and 



qualitatively explain the electrostatic mechanism of the 
current blockage. 

In Sec. IIIII we calculate the stall force F s and the stall 
effective charge q s . Our main result is that q s is pro- 
portional to the ratio of currents in blocked and opened 
channels 

q s ~ -eiVcA (4) 
Jo 

Using Eq. (Q]) we get that for ssDNA in cv-HL channel 
stall charge q s ~ — le. Intuitively this is clear because 
if the blocked current were exactly zero this would mean 
that counter ions are stuck on DNA and, therefore, com- 
pensate the DNA charge, so that the net pulling charge 
would vanish. 
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FIG. 2: The side view of a ssDNA molecule entering the chan- 
nel. 

In Sees. HVI and |V]we show that the stall charge is the 
fundamental charge so that all other charges can be ex- 
pressed in terms of it. In Sec.|lV]in order to find the cap- 
ture charge q c we calculate of the pulling force F p (X) and 
the stalling force F S {X) = —F P (X) for the partial pene- 
tration of ssDNA into the channel to the depth X < L 
(see Fig. [5]) . Then the correction to the capture barrier is 
calculated as a work which the electric field E — —V/X 
does slowly pulling DNA into the whole channel. In order 
to obtain the voltage correction —q c V to the minimum 
work necessary to overcome the capture barrier we in- 
tegrate the pulling force over X. The resulting capture 
charge q c is larger than the stall charge q s . The reason is 
that the self-energy barrier becomes smaller for a shorter 
channel. 

In Sec. |V] we discuss effective charges, which have to 
do with release rate of ssDNA, when it is trapped in the 
channel due to a hairpin in the back end (see Fig. [3]). 
There are two ways for such a DNA molecule to leave 
the channel: DNA can get unzipped by pulling electric 
field and leave to the right or DNA can escape against the 
pulling force of the electric field to the left. The former 
route dominates at large voltages, while the latter one 
dominates at smaller ones. 

It was found 0, @ that the unzipping rate exponen- 
tially grows with the voltage as 



Ru(V) = Ru(0)eM\Qu\V/k B T). 



(5) 
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We show that the "unzipping" effective charge q u — 
q s (M/N ), where M is the number of base pairs in the 
hairpin. For used in [f| M ~ 10 and Nn ~ 12, Eq. © 
gives q u ~ q s « — le in agreement with 0, Q. 

The rate of alternative escape against the voltage 
should exponentially decrease with V 



R e {V) = R e (0) exp(-q e V/k B T). 



(6) 



Here q e is the fourth effective charge, which we call 
the "escape" charge. We show in Sec. El that q e — 
\q s \(K/No), where K is number of bases in the ssDNA 
tail on the right side of the membrane when escape be- 
gins. 
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However, they cannot explain the observed dependence 
of the capture rate on the salt concentration c. In- 
deed, the persistence length of ssDNA decreases when 
c increases 16], making the conformation barrier larger, 
while in the experiment the capture rate grows with c [9( . 
So there must be another kind of barrier with the oppo- 
site c dependence. 

In Sec. I VII we suggest a mechanism for such a barrier. 
We argue that when a DNA molecule enters the channel, 
the screening cloud is squeezed in the narrow water-filled 
space surrounding the DNA. Due to this compression the 
total free energy of DNA and ions is higher for DNA in 
the channel than for DNA in the bulk. This barrier de- 
creases with c because entropy of screening atmosphere 
in the bulk solution decreases. Using even more sim- 
plified model of DNA as uniformly charge cylinder and 
the Poisson-Boltzmann approximation we show that this 
barrier is in qualitative agreement with the observed de- 
pendence R c on the salt concentration c. 

In Sec. IVIII we conclude with the summary of our re- 
sults. 



II. NEUTRALIZATION OF DNA IN THE 
CHANNEL AND THE CONTACT POTENTIAL 



FIG. 3: Unzipping of the DNA hairpin with the help of the 
voltage induced pulling force F p results in DNA translocation 
through the channel to the right. Alternatively DNA can 
escape against electric held to the left. The hairpin is shown 
schematically with the bound base pairs presented by short 
straight lines. There are M base pairs in the hairpin, Aq 
bases in the channel and K bases in the tail to the right of 
the channel. 

In Sec. IVII we are concerned with the nature of the 
capture rate barrier. Reduction of the conformation en- 
tropy due to confinement of the DNA piece in the channel 
was suggested as a natural explanation for the capture 
rate barrier Q. Indeed, ssDNA molecules in the bulk 
solution are rather flexible, with the persistence length 
about p ~ 1.4 nm 16] at 1M KC1. The channel length 
L ~ 5 nm, so it holds N p — L/p ~ 3.5 persistence lengths 
of ssDNA during translocation, and their undulations are 
restricted by the channel. The large entropic barrier due 
to this effect is NpksTAs, where As is the loss of entropy 
for one persistence length in the channel. Using As ~ 2 
we get ~ 7 ksT for this barrier. 

An additional entropy loss comes from free tails of 
DNA outside the channel [13] ■ When one end of the DNA 
chain is anchored onto the wall, the Gaussian chain has 
the free energy | In M p , where M p is the length of the tail 
in persistence lengths. The ssDNA molecules used in ex- 
periments PHH are relatively short (< 40 bases), and two 
tails can give only a barrier ~ 1— 2 fc^T. All the losses of 
conformation entropy together can explain a substantial 
part of the estimated barrier ~ 14fcsT extrapolated to 
zero voltage at salt concentration 1M KC1. 



We assume the ssDNA molecule is a rigid cylinder 
coaxial with the channel. The inner radius of the a- 
HL channel is a~0.85nm, and the radius of the ssDNA 
molecule is r ~ 0.5 nm (Fig. [1]). Salt ions are located 
in the water-filled space between them, with thickness 
6~0.35nm. The length of the channel is L~5nm. Such 
a model is even more appropriate for double helix DNA 
in a wider (say 4nm in diameter) solid state nanopore 

The dielectric constant of the channel and the ssDNA 
molecule («' ~ 2) is much smaller than that of water 
(k~80). So if ssDNA is neutralized by cations and there 
is an extra charge e at the point x located in the thin 
water-filled space between the channel internal wall, the 
electric field lines starting from this charge are squeezed 
in the thin layer (Fig. [4j. This results in a high self- 
energy U (x) of the charge [H, [H[ . 

In order to calculate U{x) we write U{x) = ecj)(x)/2, 
where <j>{x) is the electrostatic potential created by the 
extra charge at position x. If p is the distance from the 
charge e, electric field is two-dimensional at p < a (see 
Fig. and becomes uniform at larger distance p > a. 
Our numerical calculation in the limit of infinite ratio 
k/k', when all electric lines stay in the channel and at 
a/d > 1 can be well approximated by the following ex- 
pression: 



U(x) = U 1 (x)+U 2 (x) = — 
Kb 



L ( Ax 1 



1 a 

-In- 



(7) 



The origin of the two terms in Eq. {7} is illustrated in 
Fig. 0] for x — 0. At b < p < a the electric field of 
the central charge gradually spreads over all azimuthal 
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FIG. 4: An unfolded view of the water-filled space containing 
an extra K+ ion. One side of the water-filled space between 
ssDNA and the channel walls is cut by a radial half plane 
starting from the channel axis and the cut is unfolded to make 
the water filled space flat. Dashed lines represent the electric 
field lines of the charge. At p < a this electric field spreads in 
all directions and becomes uniform far from the charge. 



angles in the the whole water-filled space decaying as 
E = lej npb. This leads to the two-dimensional potential 
4>{x) = (2e//t6) hx(a/b) and produces the constant term 
U%{x) in Eq. (0. ( Eq. ([7]) works when the charge is not 
too close to the channel ends, L/2 — \x\ > a. However, 
when the charge is at channel ends Eq. |(7J) is no longer 
valid. Instead the electric field lines are attracted to the 
bulk solution, and Uq,{x) vanishes.) On the other hand, 
Ui(0) is created by the one-dimensional uniform electric 
field at distances p > a, which according to Gauss' the- 
orem is Eq — e/nab. For \x\ > the electric field at 
the closer end is stronger than that of the other end, 
therefore U\(x) decreases parabolically with |a;|, and van- 
ishes at the channel ends [III]. For L = 5 nm, a = 0.85 
nm and b = 0.35 nm Eq. ftj gives Ui(Q) = 2.9k B T and 
^2(0) = 1.7&bT where T is the room temperature. The 
total barrier U (x) of an extra K + or Cl~ ion is shown on 
Fig. [5] by the upper thick line [2(| ■ 

Now recall that there are K + ions bound to ssDNA 
phosphates in the channel. Each of them can be re- 
moved to the bulk solution creating a vacancy. The en- 
ergy penalty for this process is close to the penalty for 
placing an extra ion in the same place. Thus, energies of 
bound K + ions are —U\{x) — U% and can be shown by the 
lower full curve of Fig. [5] as a reflection of the upper one 
with respect of the x - axis. Vacancies have to overcome 
the barrier U (0) to cross the channel. Using an analogy 
with semiconductors, we can say that extra K + ions play 
the role of electrons in the conduction band, while va- 
cancies play the role of holes in the valence band. The 
peculiar result of electric field confinement in the water- 
filled space is that the energy gap 2(U\(x) + U2) has the 
maximum at x = (Fig. [5]). In the above discussion we 



ignored entropy effects II 8 1 1191 . which in principle can 
reduce self-energy barriers [2l| . 

The most important for us conclusion from above dis- 
cussion is that the large self-energy of extra charges deep 
inside the channel leads to very accurate neutralization 
of DNA by salt cations. Such nearly perfect neutraliza- 
tion was observed in computer simulations [HJ of the 
channel. 

When salt concentration in the bulk solution c is 
smaller than the characteristic concentration of K cations 
cd in DNA occupied channel, some K cations close to 
the channel ends can escape to the bulk in order to enjoy 
larger entropy in the solution. As a result there are neg- 
ative phosphate charges in the layer of width D at each 
end, and the screening (positive) charge in the adjacent 
layers of the bulk solution. These double layers (capac- 
itors) of the width D (see Fig. [1]) produce the contact 
potential —Ud, where 



Ud = k B Tln(c D /c) 



(8) 



in the channel, and prevent remaining K cations from 
leaving the channel. The contact potential appears be- 
cause negative charges in the channel are immobile (be- 
long to practically static DNA). In this sense this con- 
tact potential is similar to Donnan potential appearing 
on membranes permeable only for one sign of ions. This 
contact potential is also similar to the contact potential 
at the junction between a p-type doped and an intrinsic 
semiconductor . 

The contact potential moves down energies of both 
bands of extra cations and vacancies, while bending these 
bands up in the very ends (Fig. [5Jd). On the other hand, 
the energy band of an extra anion (shown in Fig. \Ejp by 
the dashed line) is moved up by the contact potential. 
(Without contact potential this band coincided with the 
energy band of an extra cation as shown in Fig. [5^). This 
leads to the total exclusion of anions from the channel. 
Such exclusion was also noticed in Ref. (22J. 

At c > cd both additional cations and anions freely 
enter the channel in equal number from the bulk in order 
to equilibrate their concentration. This does not lead 
to the contact potential because this process does leave 
behind layers of fixed charges. At c > cd ion current 
through the channel should be due to both extra salt 
cations and anions and, therefore, should be proportional 
to salt concentration c. On the other hand, at c < cd the 
transport is due to cations only and current is roughly 
independent on c, because Ud grows with decreasing c 
and the contact potential —Ud reduces the barrier U(0) 
(compensating for decreasing c). One can say that at 
small c transport is only due to cations neutralizing DNA. 
For example, one of them residing near the right end of 
the channel may go through channel to left end, while 
another cation from the left solution replaces it. 

These ideas provide reasonable interpretation for the 
experimental data 0] for the blocked ion Ib(c) as a func- 
tion of salt concentration c. It was found that i&(c) cx c 
at c > 1 M, while at at c < 1 M the current h(c) weakly 
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FIG. 5: Energy band diagram for K + ions (solid lines) and 
CP ions (dashed lines). The lower band represents the en- 
ergy of the cations bound to DNA phosphates (the self-energy 
necessary to create a vacancy with sign minus). The empty 
upper bands show the self-energy of the extra salt cation (solid 
line) and anions (dashed line) entering the channel, a) In the 
absence of the contact layers (c > cd). b) with contact layers 
of the width D creating contact electrostatic potential —Ud- 
Vacant phosphates are shown by empty circles. The chemical 
potential n of K + ions in the system is shown by the thin 
dotted line. 



depends on c. We interpret this data as evidence that 
Cd ~ 1 M. (See more about the experimental data for 
h(c) and their explanation in Ref. [§].) We are not trying 
here to estimate cd microscopically because this would 
require dealing with ion sizes.) 



III. EFFECTIVE CHARGE OF THE STALL 
FORCE 

The effective charge q s is defined by Eq. for the 
force F s necessary to stall ssDNA, when the ssDNA occu- 
pies the whole channel (see Fig.[T]). We show below that 
the stall charge q s is proportional to the blocked current 
lb- Let us concentrate on the case c < cd, when anions 



(CP) are excluded from the channel and all blocked ion 
current is due to cations. 

Let us assume that external electric field, E = —V/L, 
is applied in the direction opposite to x-axis. It gener- 
ates a force —NoeE on No ssDNA charges in the channel 
moving DNA along cc-axis. The opposite force N$eE acts 
on the cations. If the average drift velocity of DNA is Vd, 
and the average drift velocity of cations inside the chan- 
nel is v c we can write two momentum balance equations 
of steady state viscous motion for v c and vd- 

N eE — k c v c + k cD (v c - v D ), (9) 

-N eE + F a = k D v D + k cD (vD-v c ). (10) 

Here F a is the additional non-electrostatic force applied 
to ssDNA along — x (in the direction opposite to DNA 
motion), k c and k C D are friction coefficients of Nq cations 
with the channel walls and with DNA respectively, while 
ku is the friction coefficient of DNA with the channel 
walls. 

Although water is not included in these equations ex- 
plicitly, it is water viscosity that provides the transfer of 
momentum between moving DNA and cations and from 
them to the walls. If DNA is long (compared to Nq), the 
friction force of free ends with the water is also included 
in kr> ■ Note that voltage drops on the membrane, so that 
electric field acts only on Nq phosphates and Nq cations. 

Eq. ([9]) and (fT0|) give non-trivial predictions. First, 
when F a = 0, addition of Eq. © and (flU)) gives 

vd = -k c v c /k D . (11) 

So the average drift velocity of ssDNA is proportional 
to average drift velocity of cations. In experiment both 
v c and vd are small, so that the linear in v c and vd 
approximation is justified. Formally, Eq. (jlip looks like 
a drag of DNA by cations. Actually, this is an anti-drag, 
because the coefficient in Eq. pT|) is negative. The reason 
for the anti-drag is simple: if k c or v c were equal to zero, 
the Nq cations would transfer all momentum they receive 
from the electric field to the DNA molecule. Then DNA 
would not move at all. The transfer of negative cations 
momentum to walls makes the net force applied to DNA 
positive. Thus, DNA moves in x direction only when 
cations move in —x direction! 

Second, when F a — F s , by the definition the applied 
force stalls the ssDNA setting vd — 0. Adding Eq. © 
and p0|l at vd = we find the stalling force 

F s = k c v c , (12) 

where strictly speaking v c is the average drift velocity of 
cations at vd = 0. Actually in all experiments vd <C v c , 
which means kD is very large. Thus, one can use v c from 
experiments where F a = 0. This allows us to express v c 
through blocked current which may be written as 

h = nev Cl (13) 
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where n — Nq/L is the linear density of cations in the 
channel. Both 1^ and v c are small because of electrostatic 
barrier for ion motion in the blocked channel. Indeed, 
most of the time all cations are bound to DNA phos- 
phates. Only rarely a cation goes through middle of the 
channel contributing both into the current and the drift 
velocity v c . 

Combining Eq. (fT2"j). ^ and JT3]) we arrive at at 



-k c Ib/neE. 



(14) 



We can exclude E from Eq. (fLf|) using equation for the 
current of the open channel 



-noe 2 E/6nr]R. 



(15) 



Here = eE/6iri]R is the drift velocity of a cation in 
open channel, uq — 2cna 2 is a number of cations and 
anions per unit length of the open channel (we assume 
that they have same radius R and use Stokes formula for 
viscous resistance force because Combining Eq. 

(fH]) and HU) we get 



k c e n Q I b 
(rn-qR nig 



(16) 



The last step is to calculate k c . When the cation moves 
in the blocked channel it transfers the force eE half to 
the wall and half to the DNA. We write the force on the 
wall as 6nr]Rv c - f(R/b). If the cation size is small enough 
R << b the Stokes formula is applicable and f(R/b) = 
1/2. If the cation is large so that b — 2R << b one arrives 
at the asymptotic estimate f(R/b) = [3 In rrfg > where 
according to Ref. [23[ the coefficient (3 = 8/15. Here we 
illustrate the origin of this logarithmic expression by the 
following simple derivation leading to slightly different 
(3 = 2/3. Let us locally approximate both the DNA 
surface and internal walls of the channel as static parallel 
plane walls at z — and z — b (see Fig. ^ . Let us assume 
that center of the ball is moving in the plane z = b/2 
with velocity v c . When the ball is at x — y = and z 
axis is a polar axis of the ball let us cut the spherical 
surface of the ball by a big number of coaxial cylinders 
with the axis z and radiuses p in the range < p < R. 
We arrive at a number of rings on the surface of the ball. 
Let us estimate the force to the external wall as a sum 
over close to it rings. A distance of the ring with the 
radius p < R to the wall is h(p) = (b/2) - R + p 2 /2R, 
the local gradient of velocity is ~ v c /h, the effective area 
of the ring is 2npdp and the total momentum transfer 
rate (force) is 
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This leads to f(R/b) 



b-2R 
(17) 
2/3 for a 



/Jln^ff^ with (3 
cation in middle plane z — b/2. For a cation closer to 



2R 



external wall at finite b/2 < z < b — R we get In 

instead of In fc ff fi . Because dependence on z is only loga- 
rithmic, averaging of / over the cation positions z makes 



a negligible change to (3. Assuming that all Nq cations 
move with average drift velocity v c we find that the total 
force cations exert on the wall is 6TrNoriRv c f(R/b). This 
gives k c = 6TrNorjRf(R/b) and 



= -eN ^f(R/b). 



(18) 



For ssDNA in a-HL channel using (3 = 8/15 and a crude 
estimate In b 2 _^ R ~ 1.8 we arrive at f(R/b) ~ 1. The 
ratio no/n ~ 1 at the typical salt concentration c = 1 M. 
Thus, we arrive at Eq. (0} and q s ~ — le. 

This may look surprising because the net DNA charges 
of the DNA contact layers near the two channel ends to- 
gether are definitely larger than e. These charges, how- 
ever, do not move with DNA and do not contribute to 
the pulling DNA force and to the effective charge q s of 
the stall force. 

Up to now we dealt with the case of relatively small 
concentrations of salt, c < cp, when linear concentration 
of cations n inside the blocked channel is fixed and anions 
are excluded from the channel. One can show that at 
c > cd when current is due to both cations and anions 
in Eq. (|18p . one should drop the ratio no/n. 

We are not aware of any direct measurement of q s , 
for ssDNA in a-HL channel. In the next two sections 
we show how one can express the capture charge q c , the 
unzipping charge q u and the escape charge q e through 
the stall charge q s . 



IV. EFFECTIVE CHARGE OF THE CAPTURE 
RATE 

In the preceding section we focused on the stall force 
F s or pulling force F p = — F s in the situation, when the 
ssDNA already occupies the whole channel (Fig. [T]) . In 
order to calculate the effect of the applied voltage on the 
capture rate, we turn to the pulling force F p (X) in the 
situation, when the ssDNA penetrates only a fraction of 
the channel X/L (Fig. E2). The length of DNA in the 
channel X changes from to L while DNA enters the 
channel, and the work done by the force pulling ssDNA 



F p (X)dX, 



(19) 



This work reduces the free energy barrier of the capture 
rate (see Fig. [5]). Below we assume that all voltage drops 
on the blocked part of the channel, because unblocked 
part of the channel is so wide that its resistance is much 
smaller than that of the blocked part. Then, similarly to 
Eq. p4p . the force F p (X) can be expressed through the 
blocked current of the partially blocked state Ib(X), 



F P (X) 



kJ b (X) 



lie 



(20) 



The force F p (X) is larger than F p (L), because a shorter 
blocked channel leads to a larger current |/f,(A)| > 
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\Ib(L)\ at a given voltage V. The main reason is the 
smaller electrostatic barrier for the traversing cations. 

Combining Eq. (fr3|) , (j2"0|) and (fT4")) we arrive at a simple 
result: 

q c _ 1 f L h(X) 

In order to calculate this ratio for ssDNA in a-HL channel 
we use the experimental data [2] for the blocked current 
Ib(N) as a function of the total DNA length (in bases), 
when it is shorter than Nq — 12, namely 4 < N < 12. 
We assume that h{N) can be used for Ib{X) in Eq. (f2Tj) . 
In this way we obtain q c /q s ~ 2.5 ± 1.0 somewhat larger 
than the experimental value q c /q s ~ 1.9e/(l.le) = 1.7. 
We emphasize that in agreement with our arguments 
both theoretical and experimental values are substan- 
tially larger than e. 



and the fundamental stall charge q s . The charge q u de- 
termines the voltage dependence of the release rate of 
ssDNA, when it is trapped in the channel by an inten- 
tionally designed double helix DNA hairpin at the end 
(see Eq. ([5])). Double helix DNA of the hairpin is too 
thick to go through the a-HL channel. Thus, unzipping 
of the hairpin is necessary in order to release DNA from 
the channel to the right (Fig. [3]) . 

In experiment @ ssDNA was first inserted into the a- 
HL channel and then kept in by a relatively low voltage. 
The voltage was increased to the large value V at time 
to and probability that DNA is still in the channel at the 
time to + t was measured. This probability behaves as 
exp(— E^t) defining the rate of unzipping R u in Eq. (0. 

The rate R u was found [|| to grow with the voltage 
according to Eq. This equation defines unzipping 

charge q u . We would like to show that q u = q s (M/N ), 
where M is the number of base pairs in the hairpin. (In 
the experiment 0M=7,9, 10). 




FIG. 6: Schematic plots of several free energy barriers B(X) 
for the partial (to the length X) DNA capture. The bar- 
rier for DNA capture without electric field is nearly linear 
(dash-dotted line). Also shown are the barriers B(X) for the 
direction in which V helps to capture DNA (dashed line) and 
for the direction in which V hinders the capture (dotted line) . 
The work A(X) done by the applied potential V is shown by 
solid line. 

In Fig. [6] we illustrate the difference between capture 
barriers without a voltage, along the voltage pulling force 
and against this force. 



V. VOLTAGE ENHANCED HAIRPIN 
UNZIPPING AND ESCAPE AGAINST THE 
VOLTAGE. 

In this section we deviate from the capture rate the- 
ory and discuss relationship between unzipping and es- 
cape effective charges q u and q e defined in Introduction 




FIG. 7: Free energy W of ssDNA with a hairpin as a func- 
tion of displacement / of DNA (solid line). Positive / corre- 
spond to unzipping and release to the right, and for this case 
I is defined as the length of DNA passing by the right end 
of the channel starting from original configuration of Fig. [3] 
Negative I corresponds to the DNA escape in the direction 
of hairpin, for this case —I is defined as the length of DNA 
passing by the left end of the channel (Fig. [3}. The dotted 
line is the energy of broken base pairs, the dashed line is the 
energy W = — F p l = q s Vl/L of the voltage induced pulling 
force. Relative height of the barriers for unzipping and escape 
apparently depends on the magnitude of the voltage V. 

Indeed, unzipping would raise DNA eventually to the 
barrier Me, where e is energy of one base pair. This hap- 
pens when DNA moves along x-axis by M bases or by the 
length Imax = ML/Nq. Corresponding potential energy 
is shown in Fig. [7] by dotted line as a function of displace- 
ment I > of the DNA. The electric field reduces the 
barrier. After unzipping and moving 2M hairpin bases 
through the channel, DNA slides down the capture bar- 
rier (see Fig. [6]) and eventually reaches a constant energy 
plateau in the bulk of the right solution (Fig. [7]). 

In order to evaluate the correction to the barrier we 
should assume that DNA moves very slowly and calcu- 
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late the work of the pulling force F p = —q s (V/L), along 
the displacement l max = ML/N . This gives the correc- 



tion to the barrier F s . 
therefore 



\q s \(M/N)V = \q u \V, and, 



1u = Qs "7T" 



(22) 



M 
No' 

— le. In experiments [6j 
le. This is close to q u ~ 



We estimated above that q s 
M/Nq ~ 1 and we get q u d 
—Lie found in Ref. @. 

Our derivation is valid when the unzipping barrier cor- 
rected by voltage V is still much larger than ksT. In this 
case even with applied force, unzipping randomly alter- 
nates with zipping, while the saddle point corresponding 
to the totally unzipped hairpin is being reached in equi- 
librium way. To our mind, this condition is satisfied in 
the original experiment |6j and therefore, agreement of 
our theory with its results can be expected (24|. More 
complicated problem of unzipping of several hairpins of 
RNA was discussed recently |26j . 

Let us switch to the escape effective charge. Suppose 
DNA with a hairpin is brought to the channel by a pulling 
voltage (Fig. [3]). Let us assume that single stranded part 
of DNA has No + K bases, so that K bases have al- 
ready arrived to the bulk solution opposite to the hair- 
pin. Let us keep DNA in the channel by an applied 
voltage V, which is so small that the unzipping rate is 
smaller than the rate of escape in the direction opposite 
to the pulling force. The rate of such alternative escape 
should behave according to Eq. ([6]) and dominate at small 
enough voltages. We want to show that the escape charge 
<Ze = \ls\{K/N ). The free energy profile W for the es- 
cape process (negative I) is shown in Fig. together with 
the unzipping process (positive I). The pulling force pre- 
venting the escape is F p = \q s \V/L. DNA is climbing up 
against this force until all the DNA fits in the channel 
(no tail in the bulk solution opposite to hairpin). Start- 
ing from this point DNA slides down the capture barrier 
before reaching the energy plateau, when all DNA arrives 
to the left solution (see Fig. [5]). Therefore, displacement 



at which F p 
barrier is \l* 



works is / 



... x = —LK/Nq and the escape 
q s \(K/N )V or 



K 

^Wo- 



rn) 



The only measurement of q e we know was done in 
Ref. [H . It resulted in roughly speaking four times larger 
value of \q s \ than we estimated here, but this result 
should be taken with caution, because surprisingly in 
this experiment lb has unconventional sign. We hope 
that measurements of q e will be repeated. 



VI. ELECTROSTATIC FREE ENERGY 
BARRIER FOR DNA CAPTURE AT LOW 
VOLTAGES 



tween free energies of screening atmospheres of a rod-like 
DNA molecule in the channel and in the bulk solution, 
contributing into the capture rate barrier. We deal with 
the neutral channel, for which self-energy plays no role. 
Therefore, in the first approximation one can neglect dis- 
creteness of charges, assuming the DNA surface charges 
are uniformly smeared and the distribution of ions obeys 
Poisson-Boltzmann equation. To simplify calculation we 
look on unfolded DNA surface as an uniformly charged 
plane at z — and the unfolded inner channel wall as 
another neutral plane at z = b (Fig. We calculate 
the free energy price to push the neutral plane from in- 
finity to the distance b compressing the screening cloud. 
It is easy to verify that the price grows as the distance 
between the planes b decreases. The fact that the DNA 
surface is cylindrical is not important when the wall sep- 
aration b is smaller than the DNA radius r. We label 
the charge density of DNA surface as —a, the mean field 
electric potential as <j>. 

To find the free energy of the system we need to study 
</>(z) for a given wall separation b. Far from the channel 
ends depends only on z, and the Poisson-Boltzmann 
equation for such a system for < z < b is 



dz 2 



8iTec 



sinh (S) 



(24) 



It should be solved with the boundary conditions which 
follows from the over-all neutrality and the absence of 
charges at z < and z > b 



— \ z =b — 0, — \ Z= Q 
dz dz 



— ■ (25) 



0, k B T/e 

- 



1 - 



2 




0.5 1 1.5 2 



z , A 



FIG. 8: The potential as a function of z for b/X = 0.5, 1 and 
2 (from bottom to top), at ion concentration ce 2 /(/scr 2 ) = 1. 
Here the Gouy-Chapman length A = ksTn/ (2nae) is the unit 
of distance. 

The integration of Eq. (|24p gives 



In this section we return to the capture rate at a 
small voltage V and concentrate on the difference be- 



dz 



— cosh --- + A. 



(26) 
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With the help of the second boundary condition 
we can write A in the form A - • 1 ! 7 ,J 



Bfffc B Tc cosh (f^) , and calculate numerically 



dif> 



0(0) 



^[cosh(^)-cosh(^l)]+(^) 2 



(27) 

Numerical inversion of Eq. (|27jl gives </>(z). Then the only 
remaining parameter, the constant 0(0) in the integral 
Eq. (|27[) is determined using the first boundary condition 
(|25p . In this way we can find the potential 4>(z) for any 
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b, A 



FIG. 9: Free energy of the capture barrier per one screened 
DNA charge in the channel as a function of b (in units of 
A) for c = (solid line), ce 2 / (Ib& 2 ) = 0.1 (dotted line)and 
ce 2 /(?Bcr 2 ) = 1 (dashed line). 

given wall separation b. A few examples are shown in 
Fig. E 

We emphasize that the results obtained so far are valid 
for the neutralized main part of the ssDNA in the chan- 
nel. They are not applicable to the DNA charges near 
the end of the channel, which lose their counter ions to 
the bulk solution forming two contact layers. It is the 
contact potential —Ud that makes the potential 4>{z) in 
the neutral part of the channel negative even at z = b (see 
Fig. [5]) . This is not surprising, because through Eq. (j2"rJ|) 
the potential cj>(z) is related to c, which is in turn directly 
related to U D by Eq. ©. 

The total free energy of the system per unit area of 
DNA surface can be calculated as 



IF 



area 



^(0)+E/ 



-£-<t>{z) -k B T In — 
2 c± 



The Boltzmann distribution 



c± = cexp[- 



±e(j)(z) 
k B T 1 



c±dz. 
(28) 

(29) 



reduces Eq. ||28j) to 



W 



area 





~e(f>(z)~ 


j sinh 




[k B T\ 



e(j){z)dz - -0(0). (30) 



One can find the free energy at any given b and plot it 
choosing b = oo as the reference point (Fig. O . 

For ssDNA threading through a-HL channel with a = 
Ne/[Tr(r + a)L], the Gouy-Chapman length A - ' ' 



2-7rcre 



the salt concentration c in units of l B a 2 /e 2 , and free en- 
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ergy per area in units of k B Ta/e. Here l B 



is 



the Bjerrum length (for water at the room temperature 
l B = 0.7 ran). 



AW 



. 5 



,k B T 



n 



FIG. 10: Free energy barrier per one of screened DNA charge 
in the channel as a function of the dimensionless salt concen- 
tration n = ce 2 /(iscr 2 ), for the wall separation b = 1.4 A. 

For the system of ssDNA threading through a-HL 
channel, the Gouy-Chapman length A is about 0.25 nm, 
while the wall separation is about 0.35 nm. The free en- 
ergy barrier as a function of salt concentration in this 
example is shown by the solid line in Fig. 1101 

We see that for the salt concentration c = 1 M or 
n = ce 2 /(l B <T 2 ) = 1.37, the barrier is about AW = 
0.24k B T ■ N ~ 2Ak B T. Here N = N - 2N D is num- 
ber screened DNA phosphates in the channel and Nrj 
is number of phosphates in each contact layer. We es- 
timated that N D ~ 1 and therefore, N = 12 - 2 = 10 
in the range of salt concentrations between 0.25 M and 
c = 1 M. For smaller salt concentration c = 0.5 M and 
c = 0.25 M, the barriers are 4.3fc B T and 6.0fc B T. The 
smaller c the larger barrier. This barrier being added to 
conformation entropy barrier of ssDNA discussed in In- 
troduction may invert dependence of the total barrier on 
the salt concentration and explain the observed growth 
of the capture rate with the salt concentration. 



VII. CONCLUSION 

In this paper we evaluate effective charges of DNA re- 
sponsible for the pulling and stall forces, for the voltage 
affected capture rate, as well as for voltage induced un- 
zipping ssDNA with a hairpin and for its escape against 
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the pulling force of the electric field. The stall charge q s 
is the most fundamental one, measurements of the other 
three charges can be used to evaluate it. 

The main result of this paper is the linear equation 
connecting the stall charge with the blocked ion current 
lb- In the simplest form of Eq. (fl"8|) it is applicable only 
for relatively small concentration of salt, c < cd, when 
blocked current is due to neutralizing cations only and, 
therefore, roughly speaking is salt concentration indepen- 
dent. This equation is based only on momentum conser- 
vation, does not depend on the mechanism of the ion 
current blockage or specific model of DNA and, there- 
fore, has a high degree of universality. 

We also find a new kind of the barrier for DNA capture, 
which can explain the puzzling growth of the low voltage 
capture rate i? c (0) with the salt concentration c. We 
show that such a barrier results from squeezing of the 
screening cloud of DNA when DNA enters the channel. 

The focus of this paper is on ion channels or nanopores 



barely permitting DNA translocation. Our theory can 
be applied to ssDNA translocating through a solid state 
nanopore with diameter comparable to a-hemolysin [l2l ] . 
It should also work for a double helix DNA in a solid 
state nanopore, which diameter 2a < 3 nm only slightly 
exceeds the diameter of DNA (2 nm). For a wider pore 
with 2a > 4 nm the self-energy electrostatic barrier van- 
ishes and one can use purely hydrodynamic theory of the 
stall charge [UHzI. 
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